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ABSTRACT 



The statistics of acoustic signals propagated to long 
ranges in the ocean are investigated in detail in this 
thesis. The phase random model of multipath propagation 
is extended to include finite bandwidth and/or modulated 
sources as well as multiple source configurations. The 
theoretical analyses include the derivation of many new 
probability density functions for these new cases as 
well as for the single narrowband source. 

• 

The probability density function for A, the time 
rate of change of the level in decibels for the single 
narrowband source is derived. P: (A) depends only upon v 2 , 
the single path mean square phase rate, which can be 
related to certain ocean dynamical processes. The analysis 
of finite bandwidth and/or modulated sources reveals that 
the. amplitude and ■ amplitude rate densities (including 
P^(A)) are independent of the finite bandwidth and 
modulation effects, but the density of . the time 
rate-of-change of the multipath phase <j> is sensitive to 
these effects. Thus, fitting P: (A) to histograms from 
data to find v 2 is the preferred method for determining 
this important parameter. Bandwidth effects in can be 
neglected, however, if B << 2v where B is the signal 
bandwidth. The analysis also reveals a potentially 
powerful technique for determining parameters of the 
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modulation or bandwidth of a source from the received 
multipath signal. 

The analysis of multiple sources, applicable to noise 
problems, includes important approximations to densities 
which are intractable analytically, and would involve 
significant computer time to solve exactly. In addition 
to studies of the amplitude densities, significant progress 
has been made in solving for the amplitude rate densities 
and the joint densities of amplitude and amplitude rate. 

In addition to providing valuable confirmation of 
much of the theoretical analysis, a computer simulation 
of phase random multipath propagation also confirms that 
for M j> 4 paths phase random multipath conditions begin to 
closely approach the asymptotic conditions for M *>. 

Data at 220Hz and 406Hz received by drifting 
sonobuoys in the Atlantic at approximately 300 km in range 
were analyzed. Values of v 2 obtained support an internal 
wave model for the relevant dynamical process. The 
modulation theory uncovered a heretofore unrecognized 
modulation in the data due to an error of the Doppler 
tracking system. Predictions of crossing rates including 
this modulation effect are in good agreement with the data. 

Other data at 15Hz and 33Hz propagated to ranges 
between 250 km to 450 km in the Pacific in which deliberate 
modulation was introduced, once again provide excellent 
confirmation of the theory. Measured values of v 2 vary 
significantly from run to run and are not consistent with 
an internal wave model, indicating some other mechanism 
(i.e., tidal, rough scattering) must account for the fully 
saturated phase random nature of the data. The technique 
for determining modulation parameters was used, and for 
the 76 modulated runs analyzed, the average error in 
determining the actual bandwidth of the modulation of the 
source from the received multipath signal was 8%. 
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£ 


defined as 10 log.^ e = 4.34... 
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|p| given P = P Q (= E [p | p = p Q ] ) 
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phase crossing rate 


9 n 


single path phase for the nth path 


J v (z) 


Bessel's function of order v 


K v (z) 


the modified Bessel function of order v 


K [tt/2 , y1 


complete elliptic integral of the first kind 


L 


number of sources 


LD 


number of sources with different intensity 
and different mean square phase rates at 
the receiver 


LE 


number of sources with equal intensity and 
equal mean square phase rates at the 
receiver 


A 


level in dB's of the short time average 
mean-square pressure (A = 10 log 1Q x) 


A 


time rate of change of A 


M(t) 


temporal modulation function 


M(t) 


time rate of change of M(t) 


M (u) 
X 


characteristic function of random 
variable x 


M (to, a) 

Xi'X 2 


two-dimensional characteristic function of 
the joint pdf of random variables Xj_ 
and x 2 


M (a , 6 , z) 


the confluent hypergeometric function 
(sometimes written $(a,6,z)] 


N 


number of groups of sources, each group with 
LE i (i = 1,2 , . . . ,N) sources per group 


M 


number of independent propagation paths 
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y 


the long time average mean square pressure 
at the receiver for a single source 
(y = 2a 2 i E[ X 1) 


V 2 


the single path mean-square phase rate 
(v 2 = E [3 2 ]V n) 


p x <x) 


probability density function (pdf) of 
random variable x 


P X 1 .X 2 (X l' X 2 > 


the joint pdf of random variables X]_ 
and x 2 


p Y (x) 

N* 

<Px (S> 


the pdf of random variable x f° r N groups 

the Laplace transform of the pdf of 
random variable x 


r 


single path amplitude at the receiver 


P 


short time average root-mean-square pressure 


P 


time rate of change of p 


0 G 


variance of phase in Gaussian phase 
modulation 




one half of the long time average 
mean-square pressure at the receiver 
(c 2 = Wr 2 /2) 


W (z) 

v / y 


Whittaker's function 


X,Y 


quadrature components of the signal envelope 


X/Y 


time rate of change of the quadrature 
components 


$ 1 (a / 3 ,Y?x,y)r 


degenerate hypergeometric series in two 
variables 


X 


short time average mean-square pressure 

(x= p 2 ) 


♦ 

X 


time rate of change of x 


p 


multipath phase 



multipath phase rate 

Euler psi function [>(x) = lnF (x) J 

digamma function [^(x) = d/dx 'P (x) J 



i 

'Mx) 
4' 1 (x) 



INTRODUCTION 



Acoustic signals propagated to long ranges in the 
ocean, tens to hundreds to thousands of kilometers, via 
all modes of propagation including surface ducts, the 
deep sound channel, or sea-surface and ocean-bottom 
reflections exhibit fluctuations in amplitude and phase 
which are now recognized to be dominated by the multipath 
interference of the acoustic field. The statistics of 
these fluctuations as well as their relationships to the 
dynamics of the ocean has been one focus of recent research 
in understanding this important physical process. The 
optimum design of sonars (e.g., the receiver operating 
characteristics), underwater communications devices, and 
in fact any system which operates via acoustic 
transmission in the sea depends upon the knowledge of the 
statistical behavior of these transmissions. 

The recognition of the dominance of the multipath 
structure on the statistics, or the assertion that long 
range multipath acoustic propagation in the ocean can be 
modelled as a phase random process has been established 
only within the last ten years, although the phase random 
process or random walk problem has been under study since 
Rayleigh (1880) [1] , and is one of the classical problems 

of mathematics and physics. Bergmann (1946) [2] was among 
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the first to speculate that observed fluctuations in 
signal intensity might be attributed to the interference 
of many paths summing in random phase. Dyer (1970) [3] 

formally applied the theory of a phase random process or 
random walk problem to long range acoustic multipath 
propagation in the ocean, and was the first to investigate 
the statistics of log transformed variables. Dyer also 
showed that even in the presence of scattering randomness 
multipath interference would dominate the statistics. 

This research in fact indicated a basic shift from the 
scattering models of earlier research which are more 
appropriate for high frequencies and short ranges when 
multipath effects are less important. Dyer also proposed 
a model of distant shipping noise based upon the precepts 
of phase random multipath propagation and continued this 
research in a later paper (1973) [4] , and most recently in 

Mikhalevsky and Dyer (1978) [5] , results of the latter 

being included as part of this thesis. This model, 
appropriate to distant shipping noise, assumes the noise 
in a band is dominated by narrowband lines discrete in 
frequency. 

Mark (1972) [6] investigated the statistics of the 

multipath propagation of finite bandwidth signals. 

Employing a systems approach, he derived general expressions 
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for the mean and variance of the received signal energy 
in terms of the correlation function of the time varying 
impulse response of the medium. He also showed that the 
gamma or Erlang probability density function is often a 
good approximation to the real pdf of the received 
energy. Much earlier, Nakagami [7] had noted the 
utility of an appropriately transformed Erlang pdf in 
approximating the densities of received HF electromagnetic 
radiation undergoing rapid fading, indicating the 
broadness of scope of the phase random model and its 
general applicability. All these efforts, however, 
concentrated on the amplitude (or related quantities) of 
the signal and did not address the amplitude rate or 
phase rate of the signal. 

Longuet-Higgins (1975) [8] in connection with research 

on random sea surface waves (another phase random process) 
introduced the joint pdf's of amplitude, amplitude rate, 
phase, and phase rate as well as the marginal densities to 
the growing body of knowledge of phase random processes. 

It remained for Hamblen (1977) [9] to formally extend the 

phase random analysis to the multipath acoustic propagation 
process, incorporating the results of Longuet-Higgins [8] 
and also S.O. Rice [10] whose extensive research on noise 
statistics were also applicable to the long range acoustic 



- 26 - 



propagation problem. Hamblen [9] also established the 
dependence of the probability density functions for a 
single narrowband acoustic source on the two parameters 
a£, one-half the mean square pressure at the receiver, 
and v 2 , the single path mean square phase rate. He also 
verified the basic results with data from an ocean 
acoustic propagation experiment. 

Concurrent with the development of the phase random 
model of multipath acoustic propagation, much research 
was and is being conducted on another important aspect 
of the problem, namely to discover what ocean dynamic 
processes are the driving mechanisms and how parameters 
of these ocean dynamic processes are related to the 
parameters of ‘the acoustic field, and v 2 . Most 
notable perhaps is the recent research of Dyson, Munk, 
and Zetler (1976) [11] who have proposed a theoretical 

model relating the dynamics of internal waves in the ocean 
to the fluctuation of the acoustic field. However, this 
model appears to have serious limitations at low 
frequencies. In fact, little research has been reported 
on the low frequency cases. 

In the following paragraphs I will introduce the 
research reported on in this thesis. There are three 
basic areas in which significant progress has been made in 



